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We derive the equation of state of tracker fields, which are typical examples of freezing quintessence 
(quintessence with the equation of state approaching toward —1), taking into account of the late- 
time departure from the tracker solution due to the nonzero density parameter of dark energy fl^. 
We calculate the equation of state as a function of Q$ for constant T = VV" /{V') 2 (during matter 
era) models. The derived equation of state contains a single parameter, lU(o), which parametrizes 
f*"^ ' the equation of state during the matter-dominated epoch. We derive observational constraints on 

CNj ' W(p) and find that observational data are consistent with the cosmological constant: —1.11 < W(q\ < 

pH ! -0.96(la). 
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PACS numbers: 98.80.Cq ; 95.36.+X 



I. INTRODUCTION 

o 

! There is strong evidence that the Universe is dominated by dark energy. Moreover, the cosmological constant fits the 
current observational data quite well. However, how much is a dark energy model close to the cosmological constant? 
In order to quantify such "distance from the cosmological constant" in the dark energy theory space, we need to 
I ■ introduce a parametrization of the equation of state, w(a), which parametrizes the deviation from the cosmological 
constant, w = — 1. 

In our former study [IHll, we derived the equation of state for certain scalar field dark energy (quintessence Q/ 
k-essence [f|) models under the assumption such that the scalar field slow rolls (w ~ —1) during the matter dominated 
era. Such quintessence exhibits thawing behavior @ : the scalar field freezes during the matter era and gradually moves 
after the dark energy dominated era so that the equation of state deviates from w ~ — 1. We find our parametrization 
applies both to thawing quintessence models and to a subset of thawing k-essence models with w ~ — 1. 

However, there are quintessence models which evolve the opposite: freezing models @ with w approaching toward 
\£2 | w = — 1 so that the scalar field gradually freezes its motion. In this paper, we derive the equation of state for 
• a class of freezing quintessence models called tracker fields Q whose equation state is nearly constant during the 
matter-dominated era. For freezing models, the equation of state deviates from —1 during the matter era so that the 
, slow-roll approximation is not a good approximation. Instead, we solve the equation of motion by expanding around 
the tracker solution (the solution to which the tracker field converges from various initial conditions). The equation of 
Q\ [ state up to first order in the density parameter of dark energy was derived by @ for an inverse power-law potential 
Q. We extend the solution to higher orders in and present a useful approximation to the equation of state for 
tracker fields. Hence we now have a physically motivated parametrization of w both for thawing quintessence and for 
(a class of) freezing quintessence. 

Apart from these motivations, it is also useful to derive w(a) for more practical purposes because we do not need to 
solve the equation of motion directly for each potential. The idea has similarity, in spirit, with the slow-roll conditions: 
the existence of inflationary solutions reduces to simple conditions without having to solve the equation of motion 
directly. 

The paper is organized as follows: In Sec. 2, by perturbing the tracker equation, we derive the equation of state for 
tracker fields to all orders in fi^. In Sec. 3, we derive the observational constraints on the parameter of the equation 
of state from Type la supernovae (SNIa) data and baryon acoustic oscillations (BAO). Sec. 4 is devoted to summary. 

II. SOLVING TRACKER EQUATION 

We consider a flat universe consisting of background matter and scalar field dark energy <fi. The equation of motion 
of cf> is 

4> + 3H<j> + V'(<j>) = 0, (1) 
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where V' = dV/d(j). The equation of state w is given by 



b 2 /2 + V 

The equation of motion Eq.([T|) can be rewritten by using w Q 



where the minus sign corresponds to <j> > 0(V < 0) and the plus sign to the opposite, k 2 = 8irG, fl^(a) is the density 
parameter of dark energy, x = (1 + lo) /(l — w) and i = dhix/dhia. Tracker fields have nearly constant w initially 
and eventually evolve toward w = — 1. 

A. Tracker Solution 

Tracker fields have attractor-like solutions in the sense that a very wide range of initial conditions rapidly converge 
to a common cosmic evolutionary track: the tracker solution Q- Taking the derivative of Eq.([3]) with respect to 0, 
we obtain the so-called tracker equation [icil - fl2j 

V' 2 ~ 2(l+u;) 2(l + w)(6 + i) (l + iu)(6 + ±) (1 + w) (6 + x) 2 ' U 

where wb is the equation of state of background matter and x = d 2 \a x / d\n a 2 . Equation Q differs from the tracker 
equation in Q where the term involving f2^ is neglected which is essential in deriving the perturbation solution of w. 

Henceforth we consider the epoch after the matter-dominated era, so that we set = 0. Since during the matter- 
dominated epoch flff, is negligible and w becomes an almost constant for tracker fields, w in that epoch is written in 
terms of r as @, H 

2(r (0 ) - 1) 

W(0) = -i + 2 ( r (0) -i)' (5) 

where the zero subscript in parentheses denotes the zeroth-order solution, neglecting the contribution of dark energy 
to the expansion rate. 

B. Perturbing the Tracker Evolution 

In order to include the effect of finite f2^, we treat it as a perturbation to the zeroth-order solution and then 
extrapolate the result to the situation where 0^ is not so small when comparing the solution with the numerical 
solution. 

We define the perturbation to the zeroth-order solution fcj to be 8w, where Sw ~ S<p/<p ~ 0(Q^) from Eq. ([3]) 
Q. Keeping all terms of order in Eq. we obtain 

9 d 2 6w 5 — 6w( \ ddw 9, . , 9 . , x ^ . . . 9w 

Q ~drf~ + 2 U ~da~ + 2 (1 ~ W{Q))SW ~ 2 W(0)(1 " <4) W a ) + 9 ^ - + w (0)) dT = °- ( 6 ) 

We find that the analysis can be made simpler if T = constant so that the last term in Eq.© is vanishing. This is the 
case for inverse power-law potentials and for exponential potentials. In fact, for V = M 4 (M/<fi) a , T = (1 + a)/a, and 
also for V = M 4 e AK *, T = 1. We limit ourselves to the case when this holds so that we can solve Eq. © without using 
5<f). Note that this condition does not hold for a constant w model [lU an d for an exponential of inverse power-law 
model: V = M 4 exp(l/«0) Q. 

By approximating fi^ by the zeroth-order solution and expanding it in terms of the scale factor (or p^ {Q) / Pb) &s 
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FIG. 1: Evolution of w(a) for inverse power- law potentials: V = M 4 (M / (jj) a . The upper plot corresponds to the a = 1 case 
and the lower one corresponds to the a = 1/4 case. The solid (black) lines denote the numerical result, the dashed (blue) lines 
denote the second order solution W2(a) given by Eq. ([9]), and the dotted lines denote the first order solution w ws (a) Eq. (|10[1 . 



we find, to all orders in P0 (o) /ps or in f2^, 



w(a) = W( ) + Sw = W(o) + 
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T (l - 2i£; (0) + 4^ 0) )(1 - 3^(o) + 12^ Q) )(l - 4 W(0) + 24^ Q) )^ (a) + ' " ' 

Eq. ([8]) is our main result. 1 In the last equation, we have also expanded in terms of fi^a) using (0^(a)/(l — f^(a))) r 
^<f>( a ) n (J2m=o 0</,(a) m )™. Up to the second order in fi^, the solution becomes 



to 2 (a) = w (0) + - 



(! - w (b)) w (Q) 
- 2w (0) + Aw 2 {Q) 



(l-^ 0) )w^ 0) (8w; ( o) -1) 



(1 - 2ifl (0 ) + 4^ 0) )(1 - 3iO(o) + 12W( 0) ) 



O^(a) 



(9) 



Note that w(a) = —1 if W(o) = — 1 an d hence the cosmological constant is contained in our w(a). This w(a) (or 102(a)) 
agrees with the solution found in Q (their Eq. (33)) up to the first order in fi^,: 



w ws (a) = iO(o) 



(l-W(o)Mo) 



i — f2,A(a). 
210(0) + <W 



(10) 



J (0) 



1 The infinite series in Eq. JHJ can be written in terms of the hypergeometric functions. 
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FIG. 2: A\ 2 as a function of tu(o). 

We find the second order solution u>2 (a) Eq. (|9]) already agrees with the numerical solutions fairly well as shown in 
Fig. 1. For a < 1 (or u>( ) < —2/3), the fractional error of the equation of state between the numerical solutions and 
u>2(a), Eq. is less than 1.8%, while the error can be as large as 5% for w ws . 

III. OBSERVATIONAL CONSTRAINTS ON w (0) 

We present the observational constraints on the equation of state parameters u>(o) . We use the second order equation 
of state W2(a) Eq. ([9]) for simplicity. 2 

As observational data we consider the recent compilation of 397 SNIa, called the Constitution Set with the light 
curve fitter SALT, by Hicken et al. [l4| and the measurements of BAO from the recent SDSS data [l5[ which is now 
consistent with both earlier SDSS data [l6| and 2dF data [I?} ■ Uncertainties in the distance modulus of a supernova 
include uncertainties in light curve fitting parameters (the maximum magnitude, stretch parameter, color correction 
parameter) and due to the peculiar velocity (400 kms -1 ) as given in fl4| . 

BAO measurements from the SDSS data provide a constraint on the distance parameter A defined by 

to be A(z = 0.35) = 0.493 ± 0.017 [H. 

The x 2 curve normalized by its minimum, A% 2 = \ 2 ~ Xmim calculated from SNIa and BAO is shown in Fig. [2] 
We marginalize over f2 m to calculate the curve. The allowed range of ui(o) is narrow: —1.11 < Wm) < — 0.96(1<t), 
-1.19 < w ( o) < -0.90(2(7), -1.28 < w {a) < -0.84(3ct). 

IV. SUMMARY 

We have derived the equation of state for a class of freezing quintessence models called tracker fields whose equation 
state is nearly constant during the matter-dominated era. By solving the tracker equation perturbatively, we could 
derive a useful approximated solution to the equation of state for tracker fields (Eq. ([5])). The solutions agree with 
the numerical solutions quite accurately. 

Our solution is also useful for pragmatic purposes in that one only have to use our equation of state without solving 
the scalar field of motion numerically for any power index a. 



2 Including higher order terms does not affect the constraints. 
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Applying the solution of w(a) truncated to the second order in fi^, 102(a), to SNIa data and BAO, we find that 
the parameter ui(o) , which parameterizes the equation of state during the matter-dominated era, is constrained to lie 
near —1 and hence the cosmological constant limit of these models is consistent with the current data. 

Combining with our previous results on thawing models [IHH, we now have three parameters (wo, K, w^m) for the 
equation of state both for thawing models and for (a class of) freezing models: 



1 + w(a) 
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(K-F(a))(F(a) + l) K + (K+F(a)VF(a)-l) K 
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(thawing quintessence) , 
(freezing quintessence) , 



(12) 



where F(a) = ^Jl + (Sl^ 1 — l)a 3 . In this dark energy theory space, dark energy is close to the cosmological 

constant, which corresponds to wq = —1 irrespective of K and iDm) = — 1, to the extent such that: 3 —1.14 < wq < 
—0.92, —1.11 < W( ) < —0.96 and no constraint on K . 
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